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Constraints on modified gravity from the observed X-ray 
luminosity function of galaxy clusters 
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ABSTRACT 

We use measurements of the growth of cosmic structure, as inferred from the observed 
evolution of the X-ray luminosity function (XLF) of galaxy clusters, to constrain 
departures from General Relativity (GR) on cosmological scales. We employ the pop- 
ular growth rate parameterization, ^l^{z)'* , for which GR predicts a growth index 
7 ~ 0.55. We use observations of the cosmic microwave background (CMB), type 
la supcrnovac (SNIa), and X-ray cluster gas-mass fractions (/gas), to simultaneously 
constrain the expansion history and energy content of the Universe, as described by 
the background model parameters: fl^, w, and fJk, i.e., the mean matter density, the 
dark energy equation of state parameter, and the mean curvature, respectively. Using 
conservative allowances for systematic uncertainties, in particular for the evolution of 
the mass-luminosity scaling relation in the XLF analysis, we find 7 = 0.51^0 15 and 
fljn = 0.27 ± 0.02 (68.3 per cent confidence limits), for a flat cosmological constant 



0.44 



+0.17 
-0.15- 



0.5lt|^:i^. When in 



(ACDM) background model. Allowing w to be a free parameter, we find 7 
Relaxing the flatness prior in the ACDM model, we obtain 7 
addition to the XLF data we use the CMB data to constrain 7 through the ISW 
effect, we obtain a combined constraint of 7 = 0.45^0 :^2 the flat ACDM model. 
Our analysis provides the tightest constraints to date on the growth index. We find 
no evidence for departures from General Relativity on cosmological scales. 
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1 INTRODUCTION 

Recently, using the observed evolution of the X -ray luminos- 
ity fun ction (XLF) of massive galaxy clusters, iMantz et al.l 
l|2008l , hereafter M08) presented new constraints on dark 
energy from measurements of the growth of cosmic struc- 
ture (see also IVikhlinin et al.ll2009l ). These results are con- 
sistent with and complementary to those based on mea- 
surements of the cosmic expansion hi story, as deduced 



I2OO7I : iKomatsu et~aLl 120091 '). such measurements argue for 
a universe that is spatially flat, with most matter being cold 
and dark, and with the energy density being currently dom- 
inated by a cosmological constant, A. 

Due to severe theoretical problems associated with 
the cosmological constant, a plethora of other dark en- 
ergy models have been proposed (for a review, see 



from distances to typ e la supernovae (iRiess et al.l Il99i 




Perlmutter et all Il999l : iDavis et al.] [2OO 

aialaxy clusters (Allen et al 



iKowalski et al 



the gas mass fraction in galaxy ( 
200^ : iRapetti et aLll2005l . (200'n ). 



Frieman et al.l 20081 ) . Recent works |Carroll et al.l 



Caldwell ct al."200 7l:lHu fc Sawickill2007l : lAmin et al.1 



20061: 



20081 : 



Bcrt schingcr fc ZukinT 20081 . and^eferences therein) have 



and Bar yon Acous- 



tic Oscillations in the distr i bution of galaxies ( Cole et al] 
I2OO5I : lEisenstein et al.1 l2005l : iPercival et al] l2007f ). In com- 
bination with independ ent measurements o f the cosmic 
microwave background (|Spergel et al.l l2003l : ISpergel et al.l 
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studied the possibility that late-time cosmic acceleration 
is not driven by dark energy but rather by gravity. These 
authors propose various frameworks to test for both time 
and scale-dependent modifications to GR at late times and 
on large scales. Although, currently, modifications to GR 
may be argued to be theoretically disfavoured, it is essential 
to test for such deviations using powerful, new data that 
are now becoming available. To measure departures from 
GR on cosmological scales, experiments sensitive to the dy- 
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namical effects of GR on these scales are required. Such ex- 
periments include, for example, the cross-correlation of the 
integrated Sachs- Wolfe (ISW) effect with other matter trac- 
ers; galaxy clustering; weak gravitational lensing; or cluster 
number counts using, e.g.. X-ray selected samples. 

General Relativity has been thoroughly tested from 
laboratory to Solar System scales. However, GR has only 
just begun to be tested on cosmological scales. Several 
authors have recently investigated a simple parame- 
teri zation of the growth rate, fim(2)^ (first introduced 
by IPeeblesI llQSOl ). to test for time - dependent modifica- 



tions t o GR (see e.g. iLinder |2005l : 
Polarski fc Gannouiil booiT 



2007 



Sapone fc Amendola 



Gannouii fc Polarski 



20081: lAcauaviva et all l2008l: iBallesteros fc Riottd l20o1 



Mortonson et al.l bOOgTlThomas et al.l 120091 ). For a growth 
index, 7, of approximately 0.55, this parameterization 
accurately models the growth rate of GR. Some authors 
llNesseris fc Perivolaropoulod l2008l : iDi Porto fc Amendolal 
I2OO8I : IWeill2008l : iGonj I2OO8I ) have recently estimated con- 
straints on 7 by combining results from measurements of 
redshift space distortions and evolution in the galaxy power 
spectrum, as well as measurements of the normalization of 
the matter power spectrum, 0-8(2:), from Lyman-a forest 
data. Using current cosmic shear and galaxy clustering 
data at low redshift, iDore et al.l (|2007l ) placed constraints 
on scale-dependent modifications to GR. These authors 
constrained two phenomenological, although physically 
motivated, modifications of the Poisson equation on 
megaparsec scales (from 0.04 to lOMpc). 

In this paper, we use the XLF experiment developed by 
M08, and data from the R OSAT Brightest Cluster Sample 
rBCS: lEbeling et al.llT998l ). the ROSAT-^SO Flux Lim ited 
X-ray cluster sample fREFLEX: iBohringer et al.11200'^ '). the 
MAssive Cluster Survey (MACS: lEbeling et alTl200ll i and 
the 400 square degre es ROSAT¥SPC cluster survey (400sd; 
iBurenin et al.ll2007l ). to constrain departures from GR on 
scales of tens o f megaparsecs over th e redshif t range z < 0.9. 
We u se CMB (|Komatsu et al.ll2009l), SNIa (iKowalski et al.1 
|2008| ). and cluster /gas data (|Allen et all |2008| ) to simul- 
taneously constrain the background evolution of the Uni- 
verse. We examine three background models: fiat ACDM, 
fiat ii;CDM and non-flat ACDM. We employ a Markov Chain 
Monte Carlo (MCMC) analysis, accounting for systematic 
uncertainties in the experiments. Our results represent the 
first constraints on 7 from the observed growth of cosmic 
structure in galaxy clusters. 



2 PARAMETERIZING THE GROWTH OF 
COSMIC STRUCTURE 

In GR, the evolution of the linear matter density contrast 
5 = 5pm/ p-m, where pm is the mean comoving matter den- 
sity and 5pni a matter density fiuctuation, can be calculated 
in the synchronous gauge by solving the scale-independent 
equation 



5 + 2-5 = 4^G-Kpm5, 



(1) 



where 'dot' represents a derivative with respect to time, and 

a is the cosmi c scale factor . 

FoUowing iLahav et all ([l99j) and IWang fc Steinhardtl 

(| 19981 ) ■ several authors (see e.g. iHuterer fc Linded 20071 : 



iLinder fc £ah3[2003) have parameterized the evolution of 
the growth rate as /(a) = d\n5 / d\na = f2m(a)^. Recasting 
this expression we have the differential equation 



d5_ 

da 



S, 



(2) 



where 7 is the growth index, and f2m(a) ~ ilmCi^^ / E{a)^ . 
Here, E{a) — H{a)/Ho is the evolution parameter, H{a) 
the Hubble paramet er, and Hp its presen t-day value. It has 
been shown (see e.g. ILinder fc Cahnll2007l ) that for 7 ~ 0.55, 
Equation[2]accurately reproduces the evolution of S obtained 
from Equation[T] U sing the Einstein-Boltzmann code CAm43 
iLewis et al.ll2000l ). we find that the linear growth S{a) ob- 
tained from Equation [2] with 7 ^ 0.55 is accurate to better 
than 0.1 per cent for the relevant scales, redshifts and values 
of cosmological parameters. Therefore, we adopt 7 ~ 0.55 as 
a reference, from which to determine departures from GR. 

Equation [2] provides a phenomenological model for the 
growth of density perturbations that allows us to test depar- 
tures from GR without adopting a particular, fully covari- 
ant modified gravity theory. In the absence of such an al- 
ternative gravity theory, we perform consistency tests using 
convenient parameterizations of the background expansion, 
within GR. We investigate three expansion models that are 
well tested with current data: fiat ACDM, a constant dark 
energy equation of state, k;CD]V0, and non-fiat ACDM. We 
can write a general evolution parameter for these models 
as0 



E{a) = [: 



-I- 17k a 



1/2 



(3) 



where w = —1 for the ACDM models, Side is the cosmolog- 
ical constant/dark energy density, and fik is the curvature 
energy density, which is for fiat models. 

The growth rate Q,ni{a)'^ conveniently tends to 1 in the 
matter dominated era (high 2), thereby matching GR for 
any value of 7. Thus, we naturally match the initial value of 
5 in Equation [2] at high-2 with that of GR (see Section 



3 ANALYSIS OF THE X-RAY LUMINOSITY 
FUNCTION 

We have incorporated the growth index parameterization 
into the code developed by M08. Briefiy, in the XLF analysis, 
we compare X-ray fiux-redshift data from the cluster sam- 
ples to theoretical predictions. The relation between clus- 
ter mass and observed X-ray luminos ity is calibrated using 
deep er pointed X-ray observations (|Reiprich fc Bohringerl 
|2002| ). 



^ http://www.camb.info/ 

^ This model is only used as a expansion model, and does not 
assume the presence of dark energy. Therefore, we do not include 
dark energy density perturbations. The evolution of the density 
perturbations, due only to matter, are modelled using 7. 

Although massless neutrinos and photons are included in the 
analysis, they are neglible at late times. 
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3.1 Linear theory 



3.3 Mass— luminosity relation 



The variance of the linearly evolved density field, smoothed 
by a spherical top-hat window of comoving radius R, enclos- 
ing a mass M — AirpmR'' /S, is 



1 f°° 

o\M,z) = —j^ k^P{k,z)\WM{k)\' 



dk , 



(4) 



where Wmik) is the Fourier transform of the window func- 
tion, and P{k,z) oc k""T^{k, Zt) D{zY is the linear matter 
power spectrum as a function of the wavenumber, fc, and red- 
shift, z. Here, ris is the scalar spectral index of the primordial 
fluctuations, T{k,z%) is the matter transfer function at red- 
shift Zt, and D{z) = S{z)/5{zt) = a{M, z)/a{M, zt) is the 
growth factor of linear perturbations, normalized to unity at 
redshift Zt- We choose Zt = 30, well with in the matter dom- 
inated era (|Bertschine:er fc ZukinI [2OO8I ) . Using CAMB, we 
calculate T{k,zt) assuming that GR is valid at early times 
{z > Zt) and that modifications to GR are scale-invariant. 
For z < Zt, we calculate D{z) using Equation [2] 



3.2 Non-linear N-body simulations 



Using large N-body simulations, Ijenkins et alj (|200ll ) and 
lEvrard et aP (|2002l ) showed that the mass function of dark 
matter halos can be conveniently fitted by the expression 



M dn(M,z) , , _i 

— ,1 '/ =Aexp -lna ' + B') 



(5) 



wher e A = 0.316, B = 0.67, and e = 3.82 (|jenkins et al.1 
I2OOII ). These fit values were determined using a spherical 
overdensity group finder, at 324 times the mean matter 
density. This formula is approximately 'universal' with re- 
spect to the cosmology assumed. The universality of this 
formula has been t e sted for a wide range of cosmologies 
llKuhlen et al.1 l2005l: iLokas et al.1 12004| : iKlynin et al.1 l2003l: 



Lindor & .Jenkins' '2003'; 'Mainini ct al.' '2003'; 'Mac cio et al] 
2004: .Francis ct al. 2009; Grossi fc Sp rhigcl 2009. ). Recently 



Warren et all |2006l ) and lTinker et al] (|2008l ') have also reex 



amined this mass function using a larger suite of simulations. 

For linear scales we use the 7-model to test the growth 
rate for consistency with GR, while for non-linear scales we 
assume GEQ. As in M08, we use Equation [S] to predict the 
number density of galaxy clusters, n, at a given M and z. 



Following M08, we employ a power-law mass-luminosity 
relation, assuming self-similar evolution between the 
meiss and X-ray lumin osity, L, of massive clusters (e.g. 
iBrvan fc NormanI [l998l ') . We emphasize, however, that we 
include a generous allowance for departures from self- 
similarity, encoded in the parameter ^, 



E{z)Ma = Mo 



L 



t3 



E{z) 



{l + zf. 



(6) 



Here Ma is the cluster mass defined at an overdensity of A 
with respect to the critical density, and log Mo and /3 are 
model parameters fitted in the MCMC analysis. 

As in M08, we assume a log-normal intrinsic scatter in 
luminosity for a given mass, t), including this as a model pa- 
rameter. It is possible that this scatter may evolve with red- 
shift; we, therefore, parameterize the evolution in the scatter 
as rj{z) = ryo (1 + T?zz), where rja is the intrinsic scatter to- 
day, and 77z is a parameter that allows for linear evolution 
in the scatter. Since c urrent data are no t able to measure 
this e volution (see e.g. lO'Hara et al.ll2006l . l2007l : IChen et al.1 
[20o3), we employ the same conservative prior on r)^ used in 
M08 (see below). 

3.4 Priors and systematic allowances 

As in M08, our XLF analysis includes five parameters {A, 
C,, r]z, B, and Sbjlthat are not constrained by current data. 
We apply the same conservative priors on these parameters 
as M08. 

We apply a Gaussian prior on A, with a mean value of 
0.316, and standard deviation of 20 per cent. This devia- 
tion conservatively spans the theoretical uncertainty in the 
mass function of Equation^ For we use the uniform prior 
(—0.35, 0.35), i.e. we allow up to ~ 25 per cent change with 
respect to self-similar evolution out to z = 1. (We also in- 
vestigate results assuming strictly self-similar evolution, i.e. 
( = 0). For 7)z, we use the uniform prior (—0.3, 0.3), i.e., we 
allow up to ~ 30 per cent evolution in the scatter to z = 1. 

Departures from hydrostatic equilibrium and spheric- 
ity introduce well-known biases in mass measurements from 
X-ray data. Following M08, we assume a 25 per cent mean 
bias, and 16 per cent scatter in the bias, with 20 per cent 
systematic uncertainties in these values, as indicated by hy- 
drodynamical simulations and weak lensing data. 



^ To constrain a particular gravity model, in addition to the ex- 
pansion history and evolution of the linear density perturbations, 
we also need the mass functio n which encompasses the n on-linear 
effects of the model. Recently [Schmidt et atl | |2009^ and lSchmid^ 
( 20091') have presen ted halo mass functions for the f{R) and DGP 
( Dvali et al.ll200Cl) gravity models, respectively. The former work 
shows that for f{R) models compatible with Solar System tests, 
equation [5] provides a good fit to the mass function over the rel- 
evant mass range within the uncertainties we allow. For more 
extreme f(R) models, or for less massive clusters (for which non- 
linear chajn^leoneflfects become important) differences arise. For 
DGP, ISchmidti 120091) shows that the mass function for massive 
halos is significantly suppressed with respect to GR. Such differ- 
ences can, in principle, be constrained strongly with XLF data. 



4 INTEGRATED SACHS- WOLFE EFFECT 

Through the integrated Sachs- Wolfe (ISW) effect, the low 
multipoles of the temperature anisotropy power spectrum 
of the CMB are sensitive to the growth of cosmic structure, 
and ther efore to dark en ergy and modified gravity models 
(see e.g. iFang et al.ll2008l ). The ISW effect arises when the 
gravitational potentials of large scale structures vary with 
time, yielding a net energetic contribution to the CMB pho- 
tons crossing them. 

B and s^, are defined in Section 3.1.1 of M08. They parameterize 
the bias, and the scatter in the bias, expected for hydrostatic mass 
measurements from X-ray data. 
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Figure 1. CMB temperature anisotropy power spectra, Ci, for 
7 = 0.55 (GR; best-fit ACDM model from the five-year WMAP 
data) (solid, red line), and for values of 7 higher than GR (dashed, 
blue lines), and lower (dot-dashed, magenta lines). The lines are 
equally spaced in 7. The data are more constraining for the higher 
values of 7, which produce larger changes in the ISW effect and 
significantly worse fits to the data. Higher values of 7 monoton- 
ically increase the ISW effect while lower values first decrease it 
and then increase it (see a similar behaviour in Hu 2008). The 
circles are the binned WMAP5 data. 



On large scales, we calculate the tran sfer function due 
to the ISW effect as (|Weller fc Lewisll2003l ) 



dt e 



-T(t) 



<i)'jl [k[t - to)] 



(7) 



where 'prime' denotes a derivative with respect to the con- 
formal time t, to is the conformal time today, r is the optical 
depth to the last scattering, and jiix) is the spherical Bessel 
function for the multipole i. The total transfer function is 
the sum of Equation [7] and the transfer function from the 
last scattering surface (LSS), A,(A:) = l^f^ik) + ApsS(fc). 
From this, we calculate the anisotropy power spectrum of 
the temperature fluctuations Ci (see Figure [!}. 

For GR, we can calculate the gravitational potential (j> 
using the Poisson equation k^cj) = —iivGa^ 5pm, for which 
(j) depends only on the matter density perturbations 
Here we use the 7-model to test CMB data for consistency 
with GR, restricting ourselves to departures from GR in 
which the ISW effect can be calculated, as in GR, using the 
Poisson equation at all scalefl. 



The ISW effect is only relevant at z < 2. For this red- 
shift range, we modify CAMB to calculate the evolution of 5, 
and therefore 0', using Equation [51 and obtain 



:47rG — d(5pm[l-nm(a)^] 



(8) 



We calculate AP™(fc) at z = 2 using GR, and use these 



values as initial conditions when evolving to z 
7-model. 



5 DATA ANALYSIS 



in the 



We use the five-year WM AP CMB data (|Dunklev et all 
I2OO9I : iKomatsu et al.l |2009| . and references therein), SNIa 
data from the "Union" compil ation of iKowalski et al.l 



20081 ). the /gas measurements of lAUen et al 
XLF data f r om t he BCS {z < 0.3, 
Ebeling et all 1 19981). REFLEX (z < 0.3, 

(0.3 < 



(2008 



Bohringer etal] |2004| ). MACS 



northern 
southern 
< 0.7, 



and 
sky; 
sky; 

55 
and 



per cent sky coverage; lEbeling et al.l I2OOII . 20071) 

400sd {z < 0.9, ~ 1 per cent sky coverage; Burenin et al.l 
|20o3) cluster samples. Adopting a luminosity limit of 2.55 x 
10***/if(f ergs"^ (0.1 - 2.4keV), we use 78 clusters above 



X 10" 
-1 



from BCS, 130 



erg s ^ cm ^ 
ergs 'cm~^ from REFLEX, 34 above 
^^ergs"^cm"^ from MACS, and 30 above 0.14 x 
erg s~^ cm" 



a flux limit of 4.4 
above 3.0 x 10"^^ 
2 X 10 

10"^^ ergs~^ cm^'^ from the 400sd survey. 

We use a Metropolis MCMC algorithm to calculate 
the posterior probability distributions of our model param- 
eters. We use a modi fied version of the CoSMOMQj code 
(I Lewis fc Bridlell2002l ) that includes additional modules to 
calcu late the likelihood for the /gas experimenlQ (|Allen et al.l 
l2008h . and for the XLF experiment (M08). 

When we use fiat ACDM as background model we fit 
seven cosmological parameters: the mean physical baryon 
density, Qhh^', the mean physical dark matter density, Qch^', 
the (approximate) ratio of the sound horizon at last scat- 
tering to the angular diameter distance, 9 (which is less 
correlated with othe r parameters than Ho, as shown by 
iKosowskv et al.ll2002l ); the optical depth to reionization, r; 
the adiabatic scalar spectral index, Us', the logarithm of the 
adiabatic scalar amplitude, ln(j4s); and the growth index, 
7. For the flat wCDM background model, we have w as an 
additional parameter; and for the non-flat ACDM model, 
fik- We marginalize over seven other parameters that model 
systematic uncer tainties in the /gas analysis (see details in 
I Allen et al.ll200^ . and seven more that account for uncer- 
tainties in the XLF analysis. For the SNIa analysis, we use, 
as provided by iKowalski et al.1 (|2008l ). a covariance matrix 
that accounts for systematic uncertainties. 



^ Using CAMB we have verified this for the relevant region of pa- 
rameter space. The other two components of the Weyl tensor that 
would contribute to (f>, the anisotropic stress and ener gy flux (for 
details on these terms see IChallinor fc Lasenbviri999l) . are negli- 
gible. 

^ For a general gravity model, none of the components of the 
Weyl tensor will be negligible and, therefore modifications to the 



Poiss on equation will be required (see e.g. lHu fc Sawicki| |2007|;|Hi] 
'200^. Interestingly, for the DGP gravity model. Is^wickTetaL 
feoOTj) showed that the used of an appropiately modified Pois- 
son equation near the horizon yields an enhancement of the ISW 
effect, and thus stronger constraints with CMB data. 
* http:/ /cosmologist.info/cosmomc/ 
^ http:/ /www. stanford.edu/~drapetti/fgas_module/ 
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Figure 2. 68 and 95 per cent confidence contours in the (left panel) f2m,7 plane for the flat ACDM background model, and (right 
panel) 10,7 plane for the flat wiCDM model, from the combination of XLF, CMB, SNIa, and /gas data. These results do not include the 
(subdominant) contribution from the ISW effect (see text and Table [T] for details). 



6 RESULTS 



As shown in en et al.l (|2008l ) and references therein, the 
combination of /gas, SNIa, and CMB data places tight con- 
straints on the expansion history and The addition of 
the XLF data allows us to place tight constraints on 7. As- 
suming a flat ACDM background model, we obtain the re- 
sults shown in the left panel of Figure [2] The results for 
a fiat, constant-It) model are shown in the right panel of 
Figure [S] The marginalized values are summarized in Ta- 
ble [U For flat ACDM we find 7 = 0.5lto lt- Interestingly, 
when allowing w to be free, we obtain similar constraints 
of 7 = 0.44lg j5. We also see an anticorrelation between 
7 and w, wherein models close to a cosmological constant 
{w — —1) are most consistent with GR (7 ~ 0.55). We have 
also investigated the effect of relaxing the assumption of flat- 
ness m the ACDM model and find 7 = 0.5llo:l6, 

with no 

significant covariance between 7 and Slk- 

For massive clusters, such as those in the cur- 
rent XLF data set, self-similar evolution of the mass- 
luminosi ty relation is a well-m otivated theoretical pre- 
diction (|Brvan fc NormanI 1 19981 ). Current Chandra data 
(Mantz et al., submitted) show good agreement with this 
prediction over the redshift range spanned by the XLF data 
{z < 0.5). Under the assumption of strictly self-similar evo- 
lution (C = 0), we find 7 = O-UtHl for the fiat ACDM 
background model. 

As discussed in Section |3J the low multipoles of the 
CMB are also sensitive to the growth rate through the ISW 
effect. Even though the ISW data currently have less con- 
straining power on 7 than the XLF0, adding the ISW leads 
to ~ 17 per cent tighter constraints on 7: for the fiat ACDM 
background model we obtain 7 = 0.451q i2j ^'^^ wCDM 
7 = 0A2ta:lt' and 7 = 0.50lJJ:i5 for non-fiat ACDM. 



Note that th e cross— correlat i on of the ISW effect w ith galaxy 
surveys (see e.g. iHo et ahllioosl : iLombriser et al.ll2009l) may pro- 
vide competitive, additional constraint on 7. 



Table 1. Marginalized Icr constraints on Qm, crgi a^nd 7, using 
all the data sets without (default) or with (default -f-ISW) the 
ISW effect, for the flat ACDM (A), flat luCDM (B), and non- 
flat ACDM (C) expansion history models. For the B model we 
obtain w = —0.93 it 0.07. The results for Qm are the same with 
or without the ISW. 



(T8, 7 (default) 



0-8, 7 (default+ISW) 



A 0.27 
B 

C 0.31 



+0.02 

0.02' 
+0.02 
-0.02' 
+0.04 
0.03' 



0.82 



+0.05 
0.04' 

+0.05 



0.51 
0.44 



+0.16 
0.15 

+0.17 



,+0.04 



-+0.14 



0.76"' 



0.51^ 



0.82"' 



0.77"* 



0.42 



+0.14 
0.13 



0.50^ 



7 CONCLUSIONS 

Combining XLF, /gas, SNIa and CMB data, we have simul- 
taneously constrained the background evolution of the Uni- 
verse and the growth of matter density fiuctuations on cos- 
mological scales, allowing us to search for departures from 
CR. We parameterize the expansion history with simple 
models that include late-time cosmic acceleration, taking 
fiat ACDM as our default model, but also investigating fiat, 
constant w models, and non-fiat ACDM models. We pa- 
rameterize the growth of cosmic structure using the growth 
index, 7, which assumes the same scale dependence as GR, 
but allows for time-dependent deviations from it. 

We have performed MCMC analyses with seven (or 
eight) interesting cosmological parameters, and an addi- 
tional fourteen parameters to encompass conservative al- 
lowances for systematic uncertainties. Marginalizing over 
these allowances, we obtain the tightest constraints to date 
on the growth index. For the fiat ACDM background model. 



we measure 7 = 0.51_ 



Allowing w or Slk to be free, we 



obtain similar constraints: 7 — 0.44_ 



and 7 = 0.51 



+0.19 



respectively. Including also constraints on 7 from the ISW 
effect, we obtain 7 = 0.45^° ^ for fiat ACDM. Currently, 
we find no evidence for departures from General Relativity. 

In the near future, improved XLF analyses should pro- 
vide significantly tighter constraints on both 7 and (, al- 
lowing us to explore more complex modified gravity models 
(Rapetti et al., in preparation). Our results highlight the 
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importance of X-ray cluster data, and the potential of com- 
bined expansion history plus growth of structure studies, 
for testing dark energy and modified gravity models for the 
acceleration of the Universe. 
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